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^ , Abstract 

^ , We investigate the one-dimensional Keller-Segel model where the diffusion is replaced by a non-local operator, 

' ' ' namely the fractional diffusion with exponent < a < 2. We prove some features related to the classical two- 

I dimensional Keller-Segel system: blow-up may or may not occur depending on the initial data. More precisely 

04 ■ a singularity appears in finite time when a < 1 and the initial configuration of cells is sufficiently concentrated. 

' On the opposite, global existence holds true for a < 1 if the initial density is small enough in the sense of the 

L^/" norm. 
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"^-^ , 

. 1 Introduction 

S; 

I— I, Chemotaxis is the directed motion of cells in response to various chemical clues. It plays a key role in developmental 
biology, and more generally in self-organization of cell populations. Several categories of mathematical models have 
^ , been proposed to describe this organization process. Depending upon the level of description required, micro-, 
QQ ■ meso- or macroscopic models can be used [21] [71 [22] . Mesoscopic models consist in kinetic (scattering) equations 
, well-suited for describing the motion of bacteria such as Escherichia coli which undergo a run and tumble process 
' [SI [3] . Macroscopic models consist in parabolic (drift-diffusion) equations and are well-suited for describing motion 
, of large cells such as the slime mold amoebae Dictyostelium discoideum \10 \ fT3]. We focus on the macroscopic 
setting in this paper. 

] The so-called Keller-Segel model exhibits a very rich behaviour, emphasized by the critical mass phenomenon 

' arising in two dimensions of space. The Keller-Segel system writes in a simple formulation [16j: 

>! dtp{t,x) = Ap{t,x) -y ■ {p{t,x)Vc{t,x)) , t>0, xeM.'^ (1.1a) 

: -Acit,x)= p{t,x). (1.1b) 

■ Here p{t, x) denotes the cell density and c{t, x) denotes the concentration of the chemical attractant. The first 
contribution in the right hand side of (jl.lap expresses the tendency of cells to diffuse under their own Brownian 
motion whereas the second term expresses their tendency to aggregate due to the presence of the chemical. In two 
dimensions of space the two tendencies are evenly balanced and the global behavior of the solution depends on the 
total mass of the cells. More precisely, for M > 8tt blow-up occurs in finite time (aggregation overwhelms diffusion) 
and for M < 8tt solutions are global in time (diffusion wins the competition) [2]. 

However in one dimension of space diffusion is always stronger than aggregation and blow-up never occurs for 
systems such ([HI]) [HISO]- 
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In this paper we study the system (jl.H) in one space dimension with the cell diffusion being ruled by fractional 
diffusion. The usual Laplacian in (jl.Iap is therefore replaced by the fractional Laplacian. The non-local parabolic 
equation writes as following: 

dtp{t, x) = -h°'p{t, x) - {pd^c) , t>0, X eR (1.2a) 
-5Lc(t,x) (1.2b) 

equipped with suitable initial condition p(0, •) — po and decay conditions at infinity. For an exponent a G (0, 2], the 
positive operator A" = (— A)"/-^ is defined in Fourier variables by A"/(^) = |CI"/(0- alternative representation 
is given by: 

JyeR \x " y\ JheR 

where Ca is some normalizing factor. 

Non-local operators, and in particular the fractional Laplacian, have received a lot of attention recently [4l[5l[6]. 
In biology the motivation comes from the fact that in many cases organisms adopt Levy-flight search strategies and 
therefore dispersal is better modelled by non-local operators [U [HI [121 [ITl [18] . Focusing on the one-dimensional case 
may seem unnatural from the biological viewpoint. However we have in mind seeking a critical mass phenomenon 
as it has been derived for the two-dimensional classical Keller-Segel model (II. ip . It appears that a = c? is the critical 
fractional exponent to state such a result. Therefore it makes only sense when c? = 2 or d = 1. We partially answer 
this issue below in the latter situation. 

The system ()1.2p was first studied in [11] where it was shown that global existence holds true for 1 < a < 2 
assuming that po E H and p'^ EL?. 

We aim at providing here global existence versus blow-up results in the same spirit as for the dichotomy arising 
in the two-dimensional classical Keller-Segel system (jl.ip . More precisely, we are able to prove that solutions are 
global in time in the 'fair-competition' case a = 1, if the total mass M is assumed to be small enough. In the 
case a < 1 we show that solutions may exist globally or may blow-up depending on the initial data. We exhibit 
explicit criteria to distinguish between these two alternatives. For the case a > 1 we improve previous work [11] by 
weakening regularity hypotheses on the initial data. 

Theorem 1 (Global existence). Consider the system (|1.2p for < a < 1 with initial data po G LP°{M.) for some 
Po > 1/a. There exists a constant Ki{a) such that the condition, 

||Po||li/= < Ki{a) , 

guarantees existence of global weak solutions. 

In addition, regularizing effects act for (|1.2p . and the density belongs to any space for any positive time t > 0. 
In the case 1 < a < 2, assume po G L^''^(K) for some po > 1- Then solutions are global in time and belong to any 
U' space for all positive time t > 0. 

To complete the picture it is natural to look for blow-up results in the super-critical case. We shall prove in 
the sequel that the aggregation contribution can overcome the diffusion effect in the case a < 1 under suitable 
restrictions on the initial data. However describing the behaviour for initial data having large mass in the case 
a = 1 remains open (the constant /^2(ck) in Theorem [2] below diverges when a ^ 1). 

Theorem 2 (Blow-up). Consider the system (II. 2p for < a < 1 m one space dimension with initial data po G 
L^{{1 + \x\)dx). There exists a constant K2{a) such that the condition. 



(^J \x\pQ{x)dx^ <K2{a)Np-°' 



excludes global existence of regular solution: a singularity must appear in finite time. 

The paper is organized as follows: in Section [5] we prove global existence, beginning with a simple but not 
complete argument based on estimates. The proof is then achieved thanks to estimates inspired by [B]. 
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In Section [3] we prove blow-up of solutions. The paper is supplemented by numerical illustrations of the two 
above-mentioned phenomena. 
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2 Global existence for small initial data: proof of Theorem [T] 

We start by stating some estimate which will be widely used throughout this section. 

Proposition 2.1 (Interpolation inequality). For any exponents < a < 1 and 1 < p < +oo, the following 
Gagliardo-Nirenberg type inequality holds true: 



pP+\x)dx < C{p,a) 



Mi/c 



(2.1) 



Proof. We distinguish between the cases a < 1 and a = I. In the former we use first the Holder inequality to 
obtain: 



pP+^{x)dx < ( / pP^^^'°'\x)dx 

2 



< C{p,a) 



IIPlll/a 



where we have used the Sobolev embedding: ij"/^ i,'2/{i-a)_ 

In the case a = 1, we can use the following general result [19]: for any A, /i, s,q,r,9 G M. satisfying the following 
relations: 

^ d d d d 

i < s,q < r < oo , < < 1 , A> , ^ < , 

s r q r 



d d 



s r 



A-- + - +(1-6')U-- + - 



d d 



q r 



0. 



(2.2) 



we have, 

ll/IL.<c 11/4.,. 11/11^,^,,. 

Applying that to the particular choice: / = p^/'^, A — 1/2, ^ — Q, s — 2, r — 2{p+ l)/p, 9 — 2/r, q ~ 2/p yields the 
result. 

Observe that proceeding as above, the exponent p cannot be chosen arbitrarily (the constraint 9 > 1 forces 
p < 2). However, there is a way to extend it to any p > 1 by slightly modifying the argument: using / = p^^^, 
A = 1/2, ^ = 0, s 2, r = 2{p+l)/p, 9 = 1/ip+l), q = 2 we get: 



pP+\x)dx 



P/(2(P+1)) 



< c 



< c 



1/(P+1) 



X p/2(p+l) 

pP{x)dx] 

pP+^{x)dx 



l/2(p+l) 



Raising this inequality to the power 2(p -f 1) leads to the result. 



□ 
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2.1 A -priori I? estimates 

We complete here some existing results first derived by Escudero [IjLj . We use Gagliardo-Nirenberg type inequalities 
instead of the Sobolev inequality used in [11^ . This allows us to study a wider range of a's. We are concerned in 
this section with the global existence of the Keller-Segel system with fractional diffusion of cells when 1/2 < a < 1, 
using simple harmonic analysis estimates. This will be extended below in Section 12.21 to any < a < 1. The 
purpose of this section is to derive simply a 'priori estimates which guarantee global existence of solutions and to set 
the stage for our approach in Section Yl72\ The constraints on the exponent a here are an artefact of the method: 
in short the interpolation of L^^"- between and yields 1/2 < a. 

As it is now standard in such systems, we aim at deriving suitable U norm of the cell density. Due to the simple 
formulation of the fractional diffusion in the Fourier space variable, we opt for p = 2. We will relax this constraint 
in the next section. We have the following estimation: 



dt2 



pWIIl2= / {-h-'^p{t,x) - d^[p{t,x)d.,c{t,x))) p{t,x)dx 
Jr 

,2 , 1 



(A°'/^p{t,x)) dx + - p^{t,x)dx. 

We then apply the Gagliardo-Nirenberg inequality (Proposition [OJ for p = 2: 

[ p{t,xfdx< C{2, a) K°''^p{t)'^ ( ( p^/"{t,x)dx\ . (2.3) 
Jr \Jr J 

In the case a = 1 we obtain the decay of the norm providing that the mass is small enough: 

ilWpml. < (-7T7T^ + ^) I pHt,x)dx. (2.4) 



dt2 '"^^ '"-^ - \ C(2,1)M 2 

It follows that 11/3(^)11^2 < ||po||i2, as soon as po G L^- It is also possible to conclude without assuming pa e L^, 
by means of regularizing effects. In fact using interpolation between and it comes out that (|2.4|1 also implies 
(when the mass is small enough): 



Therefore I|p(i)||2 becomes finite in zero time. We shall come back to that later. 
In the case a < 1 the Gagliardo-Nirenberg inequality (|2.3p implies that: 

ll"'""fe-^(- C(2.a)M,)|U.. 4)//'"'-'^- 

As opposed to the case a = 1, the quantity 1|p(^)IIli/° conserved in time. Therefore we have to develop an 

alternative strategy as in [5] for the Keller-Segel in dimension d > 2, where the criterion for global existence involves 
the L'^/'^ -norm. Here we simply use the fact that L^'" can be interpolated between and if 1/2 < a < 1. As 
a consequence we have: 

Jtl ^ (^C(2,a)M--l,(t)||-- + l) //'^^'"^'^- 

Thus if the quantity M^"-! ||po||i2^" is small enough, then ||p(t)||^2 automatically decays for every time. We will 
see later that this criterion can be ameliorated, as the L^/" (before interpolation) appears to be the critical space 
for this problem (analogous to L'^^^ in the classical Keller-Segel problem). To derive this improved criterion we 
shall understand how the norms of the cell density evolve, using more refined tools for integration by parts. 

Remark 2.2. In the case 1 < a < 2, if we assume that po e H and p'q G we can work similarly as in JJJjj 
to obtain an a-priori estimate on \\Px{t)\\L2/Tg_y Then the Sobolev inequality gives a bound on 
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2.2 A priori W estimates 

Following [4l[6], the one-dimensional fractional Laplacian can be interpreted as a 'Dirichlet to Neumann problem' 
on the two-dimensional half-space (with an appropriate modification when a ^ \). Namely it is related to the 
following minimization problem. Given a function p{x) defined for x e M (and belonging to appropriate spaces, see 
[4] for details) find a function p*(a;, y) defined on M x (0, cxd) coinciding with pix) on the boundary: p*(x, 0) = p{x), 
which minimizes the weighted functional, 



J{u) 



\yu{x,y)\'y'-''dxdy. 



When a = 1 this is nothing but the harmonic extension of p to the upper half-space. The fractional Laplacian is 
then deduced from the normal derivative of p*(x, y) on the boundary {y — 0} as described below. We will strongly 
use this minimization property. 

Proposition 2.3 (Integration by parts: fractional diffusion [6]). Assume p{x) is regular, then the following estimate 
holds true: 



pP-\x)A°'p{x) dx > 



4(P - 1) 

p2 



Proof. For the sake of completeness, we recall the main lines of the proof of Proposition 12.31 We begin with the 
case a = 1 which is somewhat simpler. 



The half-Laplacian. In short, the one-dimensional half-Laplacian Ap is the normal derivative of the harmonic 
extension on the upper-half plane of p: 



where 



Ap{x) ^ -dyp^{x, 0) , 
( -Ap^{x,y) = on M x (0, cxj) , 
1 p*ix,0) = p{x) . 



Using this characterization, we are able to integrate by parts and to estimate the following diffusion contribution 
(which appears in the proof of Theorem [1] below): 



/ 



pP-^{x)Ap{x)dx = I pP-^{x,0)\/p4x,0) ■ I'dx 

^P^^^{x,y) ■ Vp^{x,y)dxdy 
|VpS^^(a;,?/)p dxdy 
> ^^fl^ / / \V{pP/^)4x,y)\'dxdy 



> 



4(P" 


1) 








1) 
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4(P- 


1) 


p2 



(2.5) 



The a/2— Laplacian. For any < a < 2 the fractional Laplacian A"p can be interpreted as follows [4]: 

AXx) - lim [-yi-"a,p,(x,2/)] , 



where 



-V • (yi-"Vp,) (a;,y) = on M x (0, cx)) 
P*{x,0) = p{x) . 
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In the same lines as (|2.5|) we are able to estimate the following diffusion contribution: 

pP-\x)A°'p{x)dx ^ ( pP-^{x,Q)y^-^\Ip,{x,Q)-vdx 



> 



I ( VpP-\x,y)-y^-^Vp,[x,y)dxdy 
Jo Jk. 

\WpPj\x,y)\Y'''dxdy 
\W{pP/\{x,y)W"'dxdy 



> 



4{p- 


1) 








1) 






4{p- 


1) 


p2 



(2.6) 
□ 



Proof of Theorem\^ The case 1 < a < 2 has already been treated in [TT] (see Remark l2.4l at the end of the proof), 
so we focus on the situation where < a < 1 in the sequel. 

L^/" is the critical space. Following the lines of 16J and [8] we estimate the evolution of the norms of the 
cell density: 



< 



pP-\t,x)A"p(t,x)dx + 

P 

, P-1 



4(p - 1) 



Ha/2 



P 



pP+\t,x) dx 
pP+\t,x)dx. 



Using Proposition 1 2 . II we obtain: 

pP+\t,x)dx<C{p,a)\\pm^,/^ (f/\t) 



2 



therefore 



dtp 



4(p - 1) P-I 

p2C(p,a)||p(i)||iVc P 



pP+\t,x)dx. 



(2.7) 



(2.8) 



Choosing in particular p = 1/a wc obtain that the L^/" norm is time-decreasing whenever UpoIIli/q is strictly 



smaller than 



4a 



C(l/Q,a) • 



Regularizing effects. We shall prove within the next lines that the cell density p{t, ■) belongs to any space 
for arbitrary positive time, provided that the initial L^" norm is finite for some pq > 1/a. The argument follows 
the main hues of [121 [HI [Sj ■ 

First we shall relax the criterion on UpoIIli/o to 



IIpoIIli/= < 



PoC{po,a) 



(2.9) 



This ensures that the L^^-norm, which is initially finite by assumption, is decreasing in time. As a consequence, we 
get the following upper-bound for any truncation A: > 0: 

||(p(t)-A:)+Lv»<|{x:p(t,:r)>fc}r-i/^"'||(p(t)-fc) + |U™ < (^-j ||po||l.o • (2.10) 
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Second, we extend the above strategy to the derivation of ||(p(t) — k{p))^\\LP for some k{p) > to be chosen 
later: 



ClL 



+C(fc,p) / {p{t,x)-k{p)Ydx + C{k,p) / {p{t,x)~k{p)Y-^ dx 



The last term can be interpolated between and . The nonlinear contribution of homogeneity p + 1 goes 
as previously, except that we shall ensure here that \\{p{t) — fc(p))+||ii/a is strictly smaller than 4/(pC(p, a)) 
independently of time. 

Introduce the notation: Yp{t) — \\{p[t) — fc(p))+||p- We have, 

dt V p C{p,a)\\{p{t) - k{p))+\\Li/c. P J Jr 

Using the following interpolation inequality: 

i-i/p 



we obtain for k{p) large enough, thanks to (|2.10p . 

^Yp{t) < -SYpitf/^P-'^ + O {Ypit)) + 0(1) , 

where 5 is a positive constant, independant of time. 

As a standard consequence, the following estimate holds true for any time t smaller than a reference time T: 

where the constant C(T) does not depend on the initial value Yp{0). These a priori estimates guarantee that the 
LP norms of p(t) {p > po) becomes finite for t > 0. □ 

Remark 2.4 (About the case 1 < a < 2). The present method can also deal with 1 < a < 2, for which global 
existence has already been proved in Ulf - In fact we shall extend accordingly Proposition \2.1\ with a/2 derivatives 
(a > 1) as following: 

2/3 „ p 



P+\x)dx< . Afi+p(i-/5), f3 = 



_ff=/2 p + a - I 

Notice that our strategy requires weaker hypotheses on the initial data ( in particular regularizing effects can be proved 
as before). 

Remark 2.5 (Intermediate asymptotics when a = 1). It is known that for the classical two-dimensional Keller-Segel 
system the cell density in space/time rescaled variables converges to a self-similar profile when mass is subcritical [31. 
The proof of this fact strongly uses the energy structure. This question is open for the one- dimensional Keller-Segel 
system with half- diffusion under consideration here. 

Recall that when only diffusion occurs (without a chemotactic coupling), such a self-similar decay holds true. This 
can be seen via the following argumentation in Fourier variables. 

First resettle time and space: u{t, y) = (1 + t)p{t, (1 + t)y), where r = log(l + t). The new equation reads: 

dru{T, y) = -Am(t, y) + dy{yu{T, y)) . 

This writes in Fourier variable as following: 

a,u(T,o = -|el"(T,C)-Ca^^i(T,0■ 
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-10 -8 -6 -4 -2 2 4 6 8 10 

space 



Figure 1: Numerical simulation of the one-dimensional Keller Segel system (|1.2p in rescaled variables with a — 1. 
The solution converges to a self-similar profile (in red). Here mass is subcritical as opposed to Fig. [2l 



Or, equivalently, 

dr {u{t, + iHr, exp(|e|)) = . 

As a consequence, 'u(t, ^) exp(|^|) can be integrated along the characteristics outgoing from 0, where u(t, 0) = M. 
This shows that u(t, ^) exp(|^|) converges to M locally in frequency. Therefore, u{T,y)/M converges to the inverse 
Fourier transform o/exp(— which is nothing hut the Cauchy density. 

Numerical simulations clearly indicate that such a statement is expected to hold true when a chemotactic con- 
tribution is added to the diffusion equation and mass is subcritical (see Fig. Qp. 



3 Blow-up: proof of Theorem [2] 

We focus in this section on the regime a < 1, for which blow-up may occur. We exhibit a criterion involving the 
mass and the first moment of the initial cell density, in the same spirit as [8] . 

Testing the fractional diffusion Keller-Segel against an adequate function cf) (regular with slow decay at infinity 
to be precised below) writes after symmetrization: 

4: I <P{x)p{t,x)dx 11 \ ('^(^) - - p{t,y))dxdy 



sgn(x - y){(j)'{x) - (j}'{y))p{t, x)p{t, y) dxdy . 



Lemma 3.1 (An auxiliary test function). Choose any < /? < 1 satisfying a-\- (3 > 1. Introduce a C°° , sub-additive, 
increasing function (f> which satisfies: (j){x) = \x\ for \x\ < 1 and <j){x) — \x\'^~^ for \x\ > 2. Denote uj^x) = —A"(j){x): 

f (f>{x + h)- (t>{x) 

= — — 

Then we have the following pointwise estimate for u: 

uj{x) KCil + lxl^-f^). 
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Proof. We split the integral into two parts: 

\^ix + h) - cj)ix)\ 



\u;{x)\ < 



< 



W<2 \h\'+'^ 



dy 



\H^ + h)-^ix)\ 



\h\>2 



dh 



i\k\<2 \h\" 
< \(t)\w^:oaC{a) + 



{x + h) + (l){x) 



\h\>2 \h\' 
|/i|i-'' + 2(/)(x 



+a 



dh 



W>2 \h\'+'^ 
< C(|(/.|vFi.- , a) + C{a, f3) + C{a)<P{x) 
<C(|<^|H.i,oo,a,/3)(l + |a:|i-^). 



dh 



□ 



Proof of Theorem\^ The proof begins with testing the Keller-Segel (|3.ip against an auxiliary function 4>\{x) = 
(j){\x)/\ where (j) satisfies the assumptions of Lemma 13. II 

4 / (i}\{x)p{t,x)dx ^ [ {~A°'(t)x{x)) p{t,x)dx [ sgn{x - y) {(t)\{x) - (j)'^{y)) p{t,x)p{t,y) dxdy . (3.2) 
™ Jm Jr * Jrxr 

Thanks to a scaling argument and Lemma l3.11 we have the following estimate: 

|A"(/.a(2:)| < A"-iC(l + 0(Ax)) . 
As a consequence we have for the first contribution in (|3.2p : 

{-A"(f»,{x))p{t,x)dx <CMA"-i+CA" / (f,x{x)p{t,x) dx . 
I Jr 

On the other hand, we can write: 







x\'-0-\x\, 



(l){x) = \x\ + R{x) , R{x) 
(jj'ix) = sgn(x) + R'{x) . 
We have clearly < C(j){x), hence: 

\R'{Xx)\ < C<j>{Xx) = CX(j)x{x) 
Therefore, we have for the second contribution in (|3.2[) : 



if \x\ < 1 
if \x\ > 2 



-i / sgn{x - y) {(l)'^{x) - (f>')^{y)) p{t,x)p{t,y) dxdy ^ / sgn{x - y) {sgn{\x) - sgn{Xy)) p{t,x)p{t,y) dxdy 

4 ./BvB 4 



< 



{{x,y):xy<0} 



sgn{x - y)R\Xx)p{x)p{y) dxdy 

p{t,x)p{t,y) dxdy + CMX I 4'\{x)p{t,x) dx . 



Observe that the symmetry assumption on the cell density p{t,x) implies the crucial point: 
/ p{t,x)p{t,y)dxdy ^ 2 ( / p{t,x)dx] ( / p{t,y)dy 

J {(x,y):xy<0} \J x<0 J \J y>0 

2 ■ 
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space 



Figure 2: Numerical simulation of the one-dimensional Keller Segel system (|1.2p with a = 1 for large mass (cell 
density is plotted in logarithmic scale). The solution clearly blows- up (final profile is plotted in red). 



We conclude the above estimates on the 'corrected' first moment Ix{t) '■— J 4>x{x)p{t,x)dx: 

^ < CMX"-^ + CX"hit) - — + CXMIxit) 
dt 4 

< ^ (CA" - \M) + C{\"+ \M) hit) . (3.3) 
4A 

We now choose A such that the terms A" and XM are well-balanced, and such that CA" — XM — —AM/2, which 
is a negative quantity. This leads to A = (/i/Af)^/(^~"\ for some constant n depending on a, /3, and the specific 
choice of the auxiliary function 0. Inequality (j3.3p rewrites: 



To finish the argumentation, let us observe that imposing a condition of the form 

^a/(i-a)j^(Q) < CMCa-")/!!-") ^ (3.4) 

yields that the quantity Ix must vanish in finite time, which is an obstruction to global existence. 

Observe finally that /a(0) < Jjg^\x\po{x)dx, hence p.4[) is satisfied if Jjg^\x\pQ{x)dx is sufficiently small. This 
completes the proof of Theorem [21 

□ 

Remark 3.2 (On the constants as a /" 1). Tracking carefully the constants in the preceding proof, it turns out that 
fj, scales like 1/(1 — a) whereas other constants are indeed of order 1. Thus criterion (j3.4p rewrites: 

/a(0)^"" < Ci-"(1 - a)"M(2-") . 

This clearly shows that the previous argument is not expected to be extended to the case a = 1. However numerical 
simulations clearly show that a critical mass is likely to occur when a — 1 (see Fig. 0). 
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